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EXECUTIVE  SUMMARY 


Human  activity  can  generate  underwater  signals  potentially  capable  of  injuring  or  altering  the 
behavior  of  mysticetes  (baleen  whales).  Little  information  exists  on  mysticete  hearing  and  this 
information  is  not  likely  to  be  obtained  in  the  near  future.  In  order  to  provide  the  Navy  with  tools 
helpful  in  predicting  the  potential  impact  of  anthropogenic  underwater  sound  on  mysticetes,  a  model 
of  humpback  whale  ( Megaptera  novaeangliae )  hearing  was  constructed.  Anatomical  indices  of 
hearing  taken  from  the  inner  ear  of  a  humpback  were  used  to  create  a  frequency-position  function 
(i.e.,  predict  the  range  of  hearing).  Auditory  sensitivity  and  frequency-position  functions  of  the  cat 
and  human  were  then  integrated  with  the  humpback  frequency-position  function  to  create  an  audio- 
gram  for  the  humpback.  The  predicted  audiogram  was  typically  mammalian  in  shape  and  suggested  a 
maximum  acoustic  sensitivity  between  2  and  6  kHz.  A  bandpass  filter  model  of  humpback  hearing, 
consisting  of  a  series  of  overlapping  pseudo-Gaussian  filters,  was  created  and  the  model  design  var¬ 
ied  via  an  evolutionary  program  to  optimize  model  sensitivity.  Agreement  between  model  sensitivity 
and  predicted  humpback  sensitivity  always  exceeded  90  percent.  The  computational  model  of  hump¬ 
back  whale  auditory  sensitivity  is  used  as  an  auditory  weighting  function  in  assessment  of  sound 
exposure. 
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INTRODUCTION 


Over  the  last  decade  the  Navy’s  use  of  high  amplitude  acoustic  signals  in  underwater  environments 
has  come  under  criticism  for  its  potential  negative  impact  upon  marine  mammals  (Henderson,  1998; 
Miller  et  al.,  2000;  Rendell  and  Gordon,  1999;  and  see  MMPA  Report  2000).  Cetaceans  (dolphins 
and  whales)  are  of  particular  concern  because  they  rely  upon  acoustic  signals  for  communication, 
foraging,  and  navigation.  Furthermore,  mysticetes  (baleen  whales)  vocalize  at  frequencies  ranging 
from  15  to  8000  Hz,  which  suggests  that  they  are  sensitive  to  frequencies  commonly  generated  by 
anthropogenic  sources  (Clark,  1990;  Herman  and  Tavolga,  1980;  Richardson  and  Wiirsig,  1997). 
Underwater  signals  generated  by  U.S.  Navy  assets  may  mask  mysticete  vocalizations,  disturb  normal 
behavior,  cause  temporary  threshold  shifts  in  hearing  (TTS),  or  cause  permanent  hearing  damage 
(Richardson  and  Wiirsig,  1997).  To  date,  little  is  known  of  the  impact  of  Navy-produced  signals  on 
mysticete  whales.  Assessment  has  been  primarily  limited  to  observations  of  responses  to  industrial 
sound  exposure  and  playback  experiments.  Unfortunately,  high  variability  within  and  between  stud¬ 
ies  has  made  reliable  conclusions  difficult  (Frankel  and  Clark,  1998;  Malme  et  al.,  1985;  Malme  et 
al.,  1988;  Maybaum,  1989;  Richardson  et  al.,  1990;  Richardson  et  al.,  1985). 

Mysticetes  are  too  large  to  maintain  in  a  controlled  environment  necessary  for  traditional  audi¬ 
ometric  tests,  but  models  of  mysticete  hearing  can  contribute  to  predictions  of  sensitivity  to  anthro¬ 
pogenic  sounds.  Houser  et  al.  (1999,  2000)  applied  evolutionary  programming  to  a  model  of  the 
dolphin  ear,  which  consisted  of  a  series  of  overlapping  frequency-domain  filters,  and  were  able  to 
match  the  model’s  sensitivity  to  the  auditory  sensitivity  of  the  dolphin.  Given  that  direct  hearing 
measurements  are  not  likely  to  be  conducted  on  any  baleen  whale  in  the  near  future,  extending  these 
types  of  computational  models  to  mysticete  audition  is  a  logical  step  towards  predicting  sensitivity  to 
U.S.  Navy  generated  sounds. 

Due  to  a  voluminous  middle  ear  cavity  and  loosely  joined  ossicles,  it  can  be  argued  that  mysticetes 
have  a  conventional  mammalian  ear  adapted  to  low  frequency  reception  (Ketten,  1997).  Methods  for 
predicting  the  frequency  range  of  hearing  from  species-specific  auditory  anatomy  (see  Greenwood, 
1990  for  review)  have  been  advanced  and  applied  to  marine  mammals  by  D.  Ketten  (Woods  Hole 
Institute,  manuscript  in  preparation).  If  the  conventional  mammalian  ear  is  assumed,  and  a  frequency- 
position  function  can  be  predicted,  psychoacoustic  and  anatomical  measures  of  hearing  from 
terrestrial  mammals  with  conventional  ears  can  be  used  to  create  a  predictive  mysticete  auditory 
threshold  function.  A  bandpass  ear-filter  model  can  then  be  constructed  and  output  optimized  to  the 
predicted  threshold  as  has  been  performed  for  the  bottlenose  dolphin  (Houser  et  al.,  1999;  Houser  et 
al.,  2000;  Roitblat  et  al.,  1993). 

This  report  describes  the  creation  of  a  bandpass  ear  model  for  the  humpback  whale  ( Megaptera 
novaeangliae),  a  medium-size  baleen  whale  with  worldwide  distribution.  Anatomical  indices  of 
hearing  derived  from  inner  ear  histology  were  used  to  generate  a  frequency-position  function. 
Sensitivity-position  information  is  lacking  for  any  mysticete.  Thus,  the  humpback  frequency-position 
data  were  integrated  with  the  known  auditory  threshold  characteristics  of  man  ( Homo  sapiens)  and 
the  domesticated  cat  ( Felis  domesticus).  The  result  is  a  predicted  audiogram  (frequency-sensitivity 
function)  for  the  humpback  whale.  Evolutionary  programming  (EP)  techniques  were  then  used  to 
optimize  the  sensitivity  of  a  bandpass-filter  model  to  the  predicted  humpback  whale  audiogram. 
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MATERIALS  AND  METHODS 


Predicting  the  Humpback  Whale  Audiogram 

An  audiometric  function  predicting  the  frequency-dependent  relative  sensitivity  of  the  humpback 
whale  was  created  on  the  assumption  that  the  humpback  ear  could  be  modeled  as  a  conventional 
mammalian  ear.  This  was  achieved  by  integrating  the  auditory  threshold  function  with  the  frequency- 
position  function  of  two  well-studied  mammals  with  conventional  ears,  the  cat  and  human,  and  map¬ 
ping  the  resulting  sensitivity-position  functions  onto  the  frequency-position  map  of  the  humpback 
whale. 

Ten  measurements  of  the  basilar  membrane  thickness  and  width  were  made  along  the  lengths  of 
basilar  membranes  obtained  from  a  humpback  whale  (D.  Ketten,  Manuscript  in  Preparation).  Ratios 
of  thickness  to  width  were  determined  (Ketten,  1993;  Ketten,  1994)  as  a  function  of  relative  position 
along  the  basilar  membrane.  Ratios  were  normalized  against  cat  thickness  to  width  ratios  and 
frequency-position  estimates  in  order  to  create  a  frequency-position  map  of  the  humpback  basilar 
membrane  (D.  Ketten,  Manuscript  in  Preparation).  A  3rd  order  exponential  was  fit  to  the  frequency 
distribution  in  order  to  generate  a  continuous  cochlear  frequency-position  function  (figure  1,  r2  =  0.99). 


Figure  1 .  Frequency  by  relative  position  along  the  length  of  the  basilar 
membrane  for  the  humpback  whale.  A  3rd  order  exponential  function 
was  fit  to  the  data  (y  =  0.08  *  exp(0.1 5x  -  (8.74e-4)x2  +  (1 ,63e-6)x3))  in 
order  to  create  the  frequency-position  function. 

Auditory  thresholds  and  cochlear  frequency-position  functions  of  the  cat  and  human  were  inte¬ 
grated  with  the  humpback  frequency-position  function  to  create  a  humpback  audiogram.  Published 
cochlear  frequency-position  functions  were  obtained  for  the  human  (Greenwood,  1990)  and  cat 
(Liberman,  1982).  For  the  human 


/(x)  =  165.4(1021jc  -1) 

and  for  the  cat 

f(x)  =  456(1021j:  -0.8) 
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where  f(x)  is  frequency  and  x  is  a  proportion  of  basilar  membrane  length.  For  all  frequencies  at  which 
hearing  has  been  tested  in  the  human  and  cat  (Fay,  1988),  the  respective  relative  position  on  the  basi¬ 
lar  membrane  at  which  that  frequency  is  predicted  to  occur  was  determined.  Relative  position  deter¬ 
minations  for  the  cat  were  limited  to  frequencies  between  100  Hz  and  60  kHz,  the  frequency  range 
covered  by  the  experimentally  determined  frequency-position  function  (Liberman,  1982). 

Threshold  intensities  (W/cm2)  of  the  cat  and  human  were  plotted  as  a  function  of  relative  basilar 
membrane  position,  converted  to  dB  re:  minimum  intensity,  and  fit  with  a  4th  order  polynomial 
(figure  2,  r2=0.66).  This  continuous  intensity-position  function  was  integrated  with  the  humpback 
frequency-position  function  to  produce  a  frequency-dependent  threshold  function  for  the  humpback. 
Thresholds  were  scaled  from  zero  to  one  in  order  to  produce  a  frequency-dependent  relative  sensitiv¬ 
ity  function. 


RELATIVE  BASILAR  MEMBRANE  POSITION 


Figure  2.  Relative  thresholds  of  the  cat  and  human  as  a  function  of 
relative  basilar  membrane  length.  The  relationship  was  determined 
by  integrating  experimentally  derived  frequency-position  functions 
for  each  species  with  their  respective  averaged  frequency-dependent 
thresholds.  The  fitted  function  (y  =  3.36  -  12.4*x  +  19.72*x2  -  19.1*x3  + 

1 1 ,75*x4)  provides  a  predictive  sensitivity-position  function. 

RESULTS 

Predicted  Humpback  Whale  Audiogram 

Frequency  by  position  along  the  basilar  membrane  was  predicted  for  ~  93%  of  the  basilar  mem¬ 
brane  length  (from  2  to  95%,  base  to  apex).  Predicted  frequency  range  was  approximately  30  Hz  to 
18  kHz.  The  predicted  humpback  whale  audiogram,  which  incorporated  an  extrapolation  of  fre¬ 
quency  by  position  to  100%  of  basilar  membrane  length,  is  shown  in  figure  3.  Threshold  is  plotted  on 
a  relative  scale,  since  we  do  not  have  an  estimate  for  absolute  sensitivity.  The  predicted  audiogram 
was  U-shaped  and  typically  mammalian  with  a  region  of  best  hearing,  defined  as  relative  threshold  < 
0.2,  ranging  from  700  Hz  to  10  kHz.  Maximum  sensitivity,  defined  as  threshold  values  <0.1,  ranged 
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from  2  to  6  kHz.  Reduction  in  sensitivity  was  ~16  dB/octave  above  10  kHz  and  ~  6  dB/  octave  below 
1  kHz.  The  most  insensitive  frequencies  occurred  at  100  Hz  and  frequencies  >  15  kHz. 


Figure  3.  The  predicted  relative  auditory  threshold  of  the  humpback 
whale.  The  threshold  function  was  determined  by  integrating  the  hump¬ 
back  frequency-position  function  (figure  1)  with  the  sensitivity-position 
function  derived  from  cat  and  human  audiometric  and  anatomic  data 
(figure  2). 


MATERIALS  AND  METHODS 


Optimized  Filter-Bank  Ear  Model  Design 

Models  of  the  humpback  ear  were  created  as  a  series  of  overlapping  bandpass  filters  with  a 
pseudo-Gaussian  shape  (Houser  et  al.,  1999)  that  delimited  the  bandpass  region  in  the  spectral  power 
domain.  Pseudo-Gaussian  filter  shapes  were  generated  with  peak  sensitivity  corresponding  to  the 
center  frequency  (p.)  of  the  filter  such  that 
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exp(-I(^)’) 

l  a 


where  x,  is  the  i,h  point  on  the  distribution  curve.  Each  filter  was  described  by  a  256  bin  vector  with 
each  bin  corresponding  to  a  100  Hz  binwidth  such  that  the  frequency  range  of  the  filter  was  0.1  to 
25.6  kHz. 


Filters  were  distributed  from  100  Hz  to  18  kHz  according  to  their  respective  center  frequency. 
Placement  of  filter  center  frequency  (jj)  was  calculated  as  a  fractional  power  of  the  frequency  range 
emulating  the  non-uniform  spacing  of  characteristic  frequencies  on  the  basilar  membrane  (Geisler 
and  Cai,  1996).  The  equation  was 

At  =  180F" 

where  Fn  is  defined  as  the  total  number  of  filters  used  in  the  model,  fj  was  the  jth  filter,  and  the  con¬ 
stant  180  was  used  to  describe  the  predicted  range  of  hearing  (i.e.,  180x  100  Hz  binwidth  =  18  kHz). 
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Filter  sensitivity  was  varied  through  an  amplitude-scaling  factor  ( S )  and  filter  width  controlled  by  a 
variable  3-dB  bandwidth  (Q3)  function.  These  factors  have  been  described  previously  (Houser  et  al., 
1999;  Houser  et  al.,  2000).  Substitution  of  parameters  into  the  pseudo-Gaussian  equation  produced 
the  final  filter  form  function: 
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Evolutionary  Programming  (EP) 

Parameters  determining  filter  shape  and  distribution  were  submitted  to  an  EP  scheme  with  self- 
adaptive  mutation  (Fogel,  1995)  and  a  Cauchy  mutation  operator  (Chellapilla  and  Fogel,  1997).  The 
EP  scheme  incorporated  a  (20+20)-EA  (evolutionary  algorithm)  such  that  populations  consisted  of 
20  “parent”  parameter  sets  with  each  “parent”  producing  1  “offspring”  parameter  set  per  generation 
(Schwefel,  1981). 

Parameters  of  the  evolutionary  program  are  given  in  Table  1  along  with  the  respective  initializa¬ 
tion  boundaries  and  standard  deviation  used  in  calculating  mutation  step-size.  Parameters,  including 
the  base  value  (y)  of  the  amplitude  scaling  factor  (S)  and  the  equation  determining  Q3,  were  mutated 
via  a  Cauchy  random  variable  (Chellapilla  and  Fogel,  1997).  The  total  number  of  filters  (Fn)  was 
mutated  in  a  probabilistic  manner  such  that  there  was  an  equal  probability  that  F„  would  increase  by 
1  or  2,  decrease  by  1  or  2,  or  stay  the  same,  if  20  <  F„  <  400.  If  Fn<  20,  there  was  an  equal  probabil¬ 
ity  that  Fn  would  increment  by  1 ,  2,  or  stay  the  same.  Conversely,  if  Fn  >  400,  there  was  an  equal 
probability  that  Fn  would  decrement  by  1,  2,  or  remain  the  same.  Thus,  minimum  and  maximum  pos¬ 
sible  values  of  Fn  were  19  and  401,  respectively. 


Table  1 .  Model  type  and  associated  parameter  values  with  initialization  limits,  initial  standard 
deviations,  and  description  of  parameter  function. 


Parameter 

Minimum 

Initialization 

Limit 

Maximum 

Initialization 

Limit 

Initial  Standard 
Deviation 

Definition  (a) 

y 

0 

10 

0.5 

base  value  for  filter  amplitude 
scaling  (used  to  calculate  S) 

m 

0 

10 

0.5 

slope  of  the  equation 
determining  Q 

b 

0 

2 

0.15 

intercept  of  the  equation 
determining  Q 

X 

0 

0.025 

0.001 

coefficient  of  the  exponent  in 
the  equation  determining  Q 

Fn 

(b) 

(b) 

(c) 

filter  number 

(a)  See  Houser  et  al.  (2000)  for  details  on  the  equations  determining  S  and  Q 

(b)  Filter  number  explicitly  set  to  40 

(c)  Probabilistic  mutation  limited  to  integer  step  sizes  of  +  2 


Following  each  generation  of  the  evolutionary  program,  defined  by  parameter  cloning  and  muta¬ 
tion,  sets  of  parameter  values  were  inserted  into  the  filter  function  to  create  a  bank  of  filters.  Each 
filter  bank  was  evaluated  for  its  sensitivity  through  a  simulated  audiometric  assessment  at  {0.1,  0.2, 
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0.3, 0.9}  and  {1.0,  2.0,  3.0, ...,  19.0}  kHz,  for  a  total  of  28  comparison  frequencies  (Houser  et 
al.,  1999).  The  filter  bank  output  formed  the  response  curve  of  the  ear  model,  which  was  normalized 
and  compared  to  the  predicted  humpback  whale  audiogram.  The  absolute  value  of  the  maximum 
deviation  between  the  filter  bank  output  curve  and  the  predicted  humpback  audiogram  curve  (MaxD) 
was  used  as  the  performance  metric  for  tournament  selection  (Goldberg  and  Deb,  1991).  Following 
sensitivity  testing  of  all  of  the  models  in  a  generation,  selection  of  parameter  sets  for  inclusion  in  the 
next  generation  was  determined  via  tournament  selection  with  a  tournament  size  of  10. 

EPs  were  run  at  the  Navy  High  Performance  Computing  Center  (SSC  San  Diego)  on  a  Hewlett- 
Packard  V2500  multi-processor  system.  The  V2500  utilized  16  440-MHz  4-way  superscalar  PA- 
8500  processors  and  16  GB  of  RAM.  Program  code  was  multithreaded  according  to  POSIX 
standards  in  order  to  take  advantage  of  the  HPC  parallel  processing  capabilities  (Norton  and 
Dipasquale,  1997;  see  Appendix  A).  Three  EP  trials  were  performed.  Trials  were  terminated  if  MaxD 
decreased  by  less  than  0.001  over  a  series  of  100  generations. 

RESULTS 

Ear  Model  Performance 

All  ear  filter  models  performed  similarly.  The  best  performing  ear  filter  model  had  a  MaxD  =  0.09 
(figure  4b),  a  slightly  better  performance  than  the  MaxD  =  0.10  achieved  by  the  other  two  models 
(figure  4a  and  4c).  Two  of  the  models  converged  upon  model  configurations  consisting  of  401  filters 
(figure  4a  and  4b)  while  the  third  utilized  263  filters  (figure  4c).  Equations  for  determining  filter  Q3 
and  the  base  value  (y)  of  the  amplitude  scaling  factor  ( S)  were,  respectively: 

(figure  4a) 

0.3  =  0.76  *  exp((1.4e-2)  *  fx)  +  0.34,  y  =  2.92 

(figure  4b) 

Qs  =  0.69  *  exp((1.5e-2)  *  fi)  +  0.23,  y  =  3.51 

(figure  4c) 

Q3  =  0.74  *  exp((1.2e-2)  *  jS)  +  0.69,  y  =  2.06 
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Figure  4.  Comparison  of  the  threshold  of  the  three  evolved  ear 
models  to  the  predicted  threshold  of  the  humpback  whale.  Triangles 
indicate  frequencies  at  which  maximum  deviations  between  the  two 
sensitivity  curves  occurred. 

Maximum  deviations  generally  occurred  above  5  kHz,  though  the  model  utilizing  a  configuration 
of  263  filters  produced  7  such  deviations  spread  across  the  predicted  range  of  hearing.  Nevertheless, 
reduction  in  sensitivity  above  10  kHz  occurred  at  a  rate  similar  to  that  predicted  by  the  target  sensi¬ 
tivity  curve.  Frequencies  of  best  sensitivity  were  at  3  and  5  kHz  for  all  models  and  there  was  a  con¬ 
sistent  sensitivity  roll-off  below  700  Hz. 

Discussion 

Humpback  whale  basilar  membrane  morphometry  has  been  combined  with  conventional  land 
mammal  psychoacoustic  and  anatomic  data  to  produce  the  first  audiometric  sensitivity  function  for 
the  humpback  whale.  The  predicted  audiogram  spans  frequencies  at  which  humpbacks  are  known  to 
vocalize — humpback  songs  span  frequencies  ranging  from  near  infrasonic  frequencies  to  over  8  kHz 
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(Helweg  et  al.,  1998;  Helweg  et  al.,  1990;  and  Payne,  1983).  Presumably  these  frequencies  lie  within 
the  range  of  humpback  hearing  thus  giving  ecological  validity  to  the  predicted  audiogram.  Though 
the  audiogram  is  plausible,  its  zone  of  best  sensitivity  (~  2  to  6  kHz)  is  somewhat  higher  than  would 
be  expected  from  the  distribution  of  frequencies  in  humpback  whale  song.  Thus,  when  implementing 
filter  models  based  upon  the  predicted  audiogram,  some  degree  of  caution  must  be  exercised  when 
interpreting  the  results. 

The  models  created  are  frequency  spectrum  filters  or  auditory  weighting  functions.  They  allow  a 
prediction  of  how  the  auditory  system  of  the  humpback  whale  attenuates  sounds  according  to  the 
sound’s  frequency-specific  components.  Thus,  it  can  be  used  to  relate  the  predicted  magnitude  of  a 
frequency  component  from  a  complex  signal  to  that  of  other  constituent  components  once  filtered  by 
the  peripheral  auditory  system.  These  models  only  make  use  of  frequency  domain  information  and 
currently  ignore  other  signal  components  such  as  time.  However,  they  serve  as  the  basis  for  the 
development  of  more  advanced  models  capable  of  incorporating  both  time  and  frequency  domains. 

Methods  of  determining  the  impact  of  anthropogenic  noise  will  likely  continue  to  rely  heavily  on 
the  use  of  playback  experiments  and  opportunistic  observations  of  mysticete  responses  to  the  expo¬ 
sure  of  man-made  sound.  Models  of  mysticete  hearing  provided  here  could  be  utilized  by  the  Navy 
to  supplement  the  results  of  such  studies  and  contextually  improve  interpretations  of  behavioral 
responses  to  sound  exposure.  Models  will  be  augmented  in  the  future  by  inclusion  of  anatomic 
information  and  biomechanical  properties  from  additional  auditory  structures  (e.g.,  middle  ear 
complex,  transfer  functions  of  the  oval  and  round  windows).  As  more  information  on  mysticete 
hearing  becomes  available,  the  effort  to  increase  the  biomimetic  nature  of  ear  models  should  continue 
until  such  time  that  absolute  auditory  thresholds  are  experimentally  determined. 
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APPENDIX  A:  SOURCE  CODE  FOR  EVOLVING  A 
HUMPBACK  WHALE  BANDPASS  EAR  FILTER  MODEL 


The  source  code  for  the  evolutionary  program  created  to  optimize  a  series  of  pseudo-Gaussian 
designed  bandpass  filters  to  the  predicted  relative  auditory  threshold  of  the  humpback  whale  may  be 
found  at  http://www.spawar.navv.mil/sti/publications/pubs/tr/1835/HumpbackEPcode.doc.  The 
following  source  code  and  header  files  are  present: 


main.C 

app.C 

app.h 

member.C 

member.h 

eval.C 

eval.h 

nrutil.C 

nrutil.h 

interp.C 

interp.h 

random.C 

random.h 

objfhs.C 

objfns.h 

The  reader  is  referred  to  Houser  et  al.  (1999)  for  procedures  on  building  the  noise  (N)  and  signal  + 
noise  (S+N)  files  used  in  testing  the  filter  models. 

Source  code  is  written  in  aCC+  and  aCC,  the  proprietary  C++  and  C  language  for  HP  UNIX  sys¬ 
tems.  When  porting  to  other  UNIX  systems,  check  the  system-specific  documentation  to  resolve 
cross-platform  incompatibilities.  Multithreading  is  achieved  through  implementation  of  the  pthread 
libraries  for  the  HP  UNIX  aCC  programming  language. 
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